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Abstract 

We calculate the self-similar longitudinal velocity correlation function, the energy spectrum and 
the corresponding other properties using the results of the Lyapunov analysis of the isotropic 
homogeneous turbulence just presented by the author in a previous work [1]. The correlation 
functions correspond to steady-state solutions of the evolution equation under the self-similarity 
hypothesis introduced by von Karman. These solutions are numerically calculated and the results 
adequately describe several properties of the isotropic turbulence. 
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I. ANALYSIS 



A recent work of the author, which deals with the Lyapunov analysis of the isotropic 
turbulence suggests a mechanism for the transferring of the kinetic energy between 
the length scales which is based on the Landau hypothesis about the bifurcations of the 
fluid kinematic equations j^. The analysis expresses the velocity fluctuation through the 
Lyapunov theory and leads to the closure of the von Karman-Howarth equation which gives 
the longitudinal velocity correlation function for two points j^, i.e. 

df Kir) fd^f Adf\ d^f, , 

+ 2Hiri + -i£ -10^^/iri(0) (1) 



dt \ dr"^ r dr J dr"^ 

where K{r), related to the triple velocity correlation function, realizes the closure of the Eq. 
([T]) through the following relation [l| 



and u is the standard deviation of the longitudinal velocity Ur, which satisfies js, 4| 



The skewness of Auj. can be expressed as ji] 

where, k{r) is the longitudinal triple velocity correlation function, related to K{r) through 

K{r) = «^ (|: + Kr) (5) 

As the result, the skewness of dur/dr is a constant which does not depend on the Reynolds 
number, whose value is H^^O) = —3/7 [1]. The other dimensionless statistical moments are 
consequentely determined, taking into account that the longitudinal velocity difference can 
be expressed as jl| 

An, ^ e + ^ ixiv' - 1) - (C' - 1)) 

Equation ([6]) arises from statistical considerations about the Navier-Stokes equations and 
expresses the internal structure of the isotropic turbulence, where ^, rj and ( are independent 



centered random variables which exhibit the gaussian distribution functions p{C), pijf) and 
whose standard deviation is equal to the unity, and ■?/' is [ll 
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The quantity R = u\t/i^ is the Taylor-scale Reynolds number, where Xt = u/ a/ {{dur / drY) 
is the Taylor-scale, whereas the function ■il){v,R) is determined as H^{r) is known. The 
parameter x is also a function of R which is given by 



;i + 2^0' (1 + x^)) 



3/2 7 



(8) 



with ijjQ = iIj{R,0) and ipo = 1-075 1|. From Eqs. (EI) and ([71), all the absolute values of 
the dimensionless moments of Aur of order greater than 3 rise with R, indicating that the 
intermittency increases with the Reynolds number. 

The PDF of Aur can be formally expressed with the Frobenious- Perron equation 



F{Au'r)= [ [ [p{Op{v)p{C)H^Ur-Au'r)d^dr]dC 



(9) 



where 5 is the Dirac delta, whereas the spectrums E{k) and T(k) are calculated as the 
Fourier Transforms of / and -ftT jd], respectively, i.e. 
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II. SELF-SIMILARITY 

In this section the properties of the self-similar solutions of the von Karman-Howarth 
equation are studied. 

Far from the initial condition, it is reasonable that the mechanism of the cascade of 
energy and the effects of the viscosity act keeping / and E{k,) similar in the time. This is 
the idea of self-preserving correlation function and turbulence spectrum which was originally 
introduced by von Karman (see ref. [sl and reference therein). 
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In order to analyse this self-similarity, it is convenient to express / in terms of the 
dimensionless variables f = r/Ay and i = tu/Xx, i.e., / = f{i,r). As the result, Eq. 
reads as follows 



dfXTdftu\ df fdXr W^l^V H/ (11) 



di u dt \Xt J df u dt V 2 df R\ f df J R' 
where c}^//9f^(0) = — 1. This is a non-linear partial differential equation whose coefficients 
vary in time according to the rate of kinetic energy 



dt Xl 



(12) 



If the the self-similarity is assumed, all the coefficients of Eq. ffTTl) must not vary with the 

nn 

time [3|, 15|, thus one obtains 

R = const (13) 

Xt d f tu\ . 
«i = — \ ^]= const (14) 



u dt \Xt 



02 = — T~ = const, (15) 
u dt 



As the consequence of Eqs. ( 1121) and (]T3l) . Xt and u will depend upon the time according to 
XT{t) = XT{0)Jl + lOu/XUO)t, «W = ^====F=- (16) 

A/i + iUZ^/ Xrp([j ) t 



From these expressions, the corresponding values of ai and 02 are 

_ Xt d / tu\ 1 



u dtyXTj 1 + lOz/t/A: 



2 ' 



(17) 
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02 = -- 



u dt R 

The coefficient ai decreases with the time and for t — * 00, Oi — 0, whereas 02 remains 
constant. Therefore, for t — > 00, one obtains the self-similar correlation function /(f), which 
does not depend on the initial condition and that obeys to the following non-linear ordinary 
differential equation 



i^^: /!-/ df 2 fd'f 4df\ 10 .^g. 
R df^ \ 2 df RKdf^ fdf R^ 



The first term of Eq. (fT8|) represents tlie variations in time of / which does not influence 
the mechanism of energy cascade. This term is neghgible with respect to the second one 
only if f << i?, and is responsible for the asymptotic behavior of / which is expressed by 
/ ^ This behavior determines that all the integral scales of / diverge and that / 

does not admit Fourier transform, thus the corresponding energy spectrum is not defined. 
According to von Karman , we search the self-similar solutions over the whole range of 
f, with the exception of the dimensionless distances whose order magnitude exceed R. This 
corresponds to assume the self-similarity for all the frequencies of the energy spectrum, but 



for the lowest ones 



J, Isj. As the result, the first term of Eq. ( fT8|) can be neglected with 



respect to the second one, and the equation for / reads as follows 



f df 2 /dV 4^y\ iQ 



2 df R \ df"^ fdf) R' 

The analysis of Eq. ([19]) shows that / ~ 1 - + (10 + i?)/112 in the vicinity of the 
origin and that / — 1 ~ f"^^'^ when the first term of Eq. fll9p is about constant, whereas for 
large f, / exponentially decreases. Thus, all the integral scales of / are finite quantities and 
the energy spectrum is a definite quantity whose integral over the Fourier space gives the 
turbulent kinetic energy. 



III. RESULTS AND DISCUSSION 

In this section we calculate the solutions of Eq. f|T9|) and study the corresponding prop- 
erties of the self-similar solutions. To determine / and its energy spectrum, consider the 
following initial condition problem with respect to the dimensionless separation distance f 

— -F 
df 

(20) 



df \^2 V 2 fj 

This ordinary differential system arises from Eq. f[T^ and its initial condition is /(O) = 1, 
F(0) = 0. 

Several numerical solutions of Eqs. (120!) were calculated for different Taylor scale Reynolds 
numbers by means of the fourth-order Runge-Kutta scheme of integration. The cases here 
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FIG. 1: Longitudinal correlation function for different Taylor-Scale Reynolds numbers. 



analyzed correspond to i?=100, 200, 300, 400, 500 and 600. The fixed step size of the 
integrator scheme is selected on the basis of the asymptotic stability condition Af = \/2/R 
Gj, which also provides a fairly accurate description of the energy spectrum at the large 
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FIG. 2: Longitudinal triple correlation function for different Taylor-Scale Reynolds numbers. 
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wave-numbers. 

Figures [1] and [2] show the numerical solutions of Eqs. fl20l) . where double and triple 
longitudinal correlation functions are represented in terms of f, for the different values of 
R. Due to the mechanism of energy cascade, the tail of / rises with R and the maximum of 
\k\ gives the entity of this mechanism. This value is slightly less than 0.05 and agrees quite 
well with the numerous data of the literature which concern the evolution of the correlation 
functions. It is apparent that the spatial variations of k correspond to dimensionless scales 
f whose size increases with R. 

Figures [3] and H] show the plots of E{k,) and T(k) for the same Reynolds numbers. As 
the consequence of the mathematical properties of /, the energy spectrum behaves like 
E{k,) = 0{k,^) in proximity of the origin, and after a maximum is about parallel to the 
—5/3 Kolmogorov law (dashed line in Fig. [3]) in a given interval of the wave-numbers. This 
interval defines the inertial range of Kolmogorov, and its size increases with R. For higher 
wave-numbers the energy spectrum rapidly decreases with a slope which depends on the 
behavior of / in proximity of the origin and thus on the Reynolds number. 

Since K does not modify the kinetic energy of the flow, according to Eq. ([2]), the integral 
of T(/t) over the Fourier wave-numbers results to be identically equal to zero at all the 
Reynolds numbers. 
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FIG. 3: Turbulent Energy Spectrum for different Taylor-Scale Reynolds numbers. 
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FIG. 4: "Transfer function T{k)" for several Taylor-Scale Reynolds numbers. 

In the Figs. [5] and [6], skewness and flatness of Aur are shown in terms of f for the same 
values of R. The skewness H^, is flrst calculated according to Eq. (jlj) and thereafter the 
flatness if 4 has been determined using Eq. For a given R, {H^l starts from 3/7 at the 
origin, then decreases to small values, while starts from values quite greater than 3 at 
r = 0, then reaches the value of 3 (faster than tending to zero). Although H^lO) does 
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FIG. 5: Skewness of Aur at different Taylor-Scale Reynolds numbers. 
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FIG. 6: Flatness of Aur at different Taylor-Scale Reynolds numbers. 

not depend upon R, H^{f) is a rising function of R and, in any case, the intermittency of 
/S.Ur increases with R according to Eqs. (I6l) and ([7]). 

Next, the Kolmogorov function Q{r) and Kolmogorov constant C, are determined using 
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FIG. 7: Kolmogorov function for several Taylor-Scale Reynolds numbers. 
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R C 

100 1.8860 
200 1.9451 
300 1.9704 
400 1.9847 
500 1.9940 
600 2.0005 

TABLE I: Kolmogorov constant for different Taylor-Scale Reynolds number, 
the previous results. According to the theory, the Kolmogorov function, defined as 

«M^-M pi) 

is constant with respect to r, and is equal to 4/5 as long as r/Ay = 0(1). As shown in Fig. 
[71 Q{r) exhibits a maximum for f = 0(1) and quite small variations for higher f, as the 
Reynolds number increases. This maximum increases with i?, and seems to tend toward the 
limit 4/5 prescribed by the Kolmogorov theory. 

The Kolmogorov constant O, defined by E{k) Ce^^"^ / k^^^ , is here calculated as 

£;(fi:)/t5/3 

Ke(o,oo) 

where e = —3/2 du^/ dt is the rate of the energy of dissipation. In the tablelH the Kolmogorov 
constant is reported in terms of the Taylor-scale Reynolds number. The obtained values of 
C and Qmax are in good agreement with the corresponding values known from the various 
literature. 

The spatial structure of Am^, expressed by Eq. ([6]), is also studied with the previous 
results. 

According to the various works j^, [o], 10], Am^ behaves quite similarly to a multifractal 
system, where Au^ obeys to a law of the kind Aur{r) ~ r'^ in which g is a fiuctuating ex- 
ponent. This implies that the statistical moments of Aur(r) are expressed through different 
scaling exponents C(^) whose values depend on the moment order n, i.e. 

((AM,)"(r)) = Ar^(") (23) 

In order to calculate these exponents, the statistical moments of Am^ are first calculated 
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C= max ^^^^5^ (22) 
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FIG. 8: Statistical moments of Au.,. in terms of the separation distance, for i?=600. 

using Eqs. (El) and (I9j) for several separation distances. Figure M shows the evolution of the 
statistical moments of Aur in terms of f, in the case of R = 600. The scaling exponents 
of Eq. (!23|) are identified through a best fitting procedure, in the intervals (ri,f2), where 
the endpoints fi and f2 have to be determined. The calculation of ({n) and An is carried 
out through a minimum square method which, for each moment order, is applied to the 
following optimization problem 

Uan),Ar,)= r(((AM,r)-A„r^("))2rfr = min, n = l,2,... (24) 
J fi 
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C(ii) 


C(12) 
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100 


0.35 


0.70 


1.00 


1.30 


1.56 


1.82 


2.06 


2.31 


2.53 


2.76 


2.97 


3.18 


3.39 


3.59 


3.79 


200 


0.35 


0.71 


1.00 


1.29 


1.55 


1.81 


2.05 


2.28 


2.50 


2.72 


2.93 


3.14 


3.33 


3.53 


3.73 


300 


0.35 


0.71 


1.00 


1.29 


1.55 


1.81 


2.05 


2.28 


2.50 


2.73 


2.93 


3.14 


3.34 


3.54 


3.73 


400 


0.35 


0.71 


1.00 


1.29 


1.55 


1.81 


2.04 


2.28 


2.50 


2.72 


2.93 


3.13 


3.33 


3.53 


3.72 


500 


0.35 


0.71 


1.00 


1.29 


1.55 


1.81 


2.04 


2.28 


2.50 


2.72 


2.93 


3.13 


3.33 


3.53 


3.73 


600 


0.35 


0.71 


1.00 


1.29 


1.55 


1.81 


2.05 


2.28 


2.51 


2.73 


2.94 


3.15 


3.35 


3.55 


3.75 



TABLE II: Scaling exponents of the longitudinal velocity difference for several Taylor-Scale 
Reynolds number. 
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FIG. 9: Scaling exponents of Auj. for several R. Solid symbols are for the present data. Dashed 
line is for Kolmogorov K41 data Ml. Dotted line is for Kolmogorov K62 data Continuous line 

Q 

is for She-Leveque data [9( 

where (((AUf.)")) are calculated with Eqs. ([6]), fi is assumed to be equal to 0.1, whereas 
^2 is taken in such a way that ({3) = 1. The so obtained scaling exponents are shown in 
Table ([Tll) in terms of the Taylor scale Reynolds number, whereas in Fig. [9 (solid symbols) 
these exponents are compared with those of the Kolmogorov theories K41 3 (dashed line) 
and K62 (dotted line), and with the exponents calculated by She-Leveque sl (continuous 
curve). Near the origin ({n) ^ n/S, and in general the values of ({n) are in good agreement 
with the She-Leveque results. In particular the scaling exponents here calculated are lightly 
greater than those by She-Leveque for n > 8. 

The PDFs of dur/df are determined by means of Eqs. ([9]) and ([6]). Specifically, the PDF 
is calculated with direct simulations, where the sequences of the variables ^, t] and C, are 
first determined by a gaussian random numbers generator. The distribution function is then 
calculated through the statistical elaboration of the data obtained with Eq. ([6]). The results 
are shown in Fig. [TOb and [TOb in terms of the dimensionless abscissa 

dur/df 
{{dur/drfyl^ 

These distribution functions are normalized, in order that their standard deviations are 
equal to the unity. The figure represents the PDF for the several i?, where the dashed curve 
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FIG. 10: Probability distribution functions of the longitudinal velocity derivative for the different 
Taylor-Scale Reynolds numbers 

represents the gaussian distribution functions. In particular, Fig. [TOb shows the enlarged 
region of Fig. [TOk . where 5 < s < 8. According to Eq. ([6]), the tails of PDFs change with 
R in such a way that the intermittency of dur/df rises with the Reynolds number. 



IV. CONCLUSIONS 



The obtained self-similar solutions of the von Karman-Howarth equation with the pro- 
posed closure, and the corresponding characteristics of the turbulent flow are shown to be 
in very good agreement with the various properties of the turbulence from several points of 
view. 

In particular: 

• The energy spectrum follows the Kolmogorov law in a range of wave-numbers whose 
size increases with the Reynolds number. 

• The Kolmogorov function exhibits a maximum and relatively small variations in prox- 
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imity of r = 0{Xt)- This maximum value rises with the Reynolds number and seems 
to tend toward the limit 4/5, prescribed by the Kolmogorov theory. 

• The Kolmogorov constant moderately varies with the Reynolds number with an aver- 
age value around to 1.95 when R varies from 100 to 600. 

• The scaling exponents of the moments of velocity difference are calculated through a 
best fitting procedure in an opportune range of the separation distance. The values of 
these exponents are in good agreement with the results known from the literature. 

• The intermittency of the longitudinal velocity difference rises with the Reynolds num- 
ber. 

These results represent a further test of the analysis presented in Ref. ^ which adequately 
describes many of the properties of the isotropic turbulence. 
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